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Abstract. The mass loss rates, expansion velocities and dust-to-gas density ratios from millimetric observations of 1 19 carbon- 
rich giants are compared, as functions of stellar parameters, to the predictions of recent hydrodynamical models. Distances 
and luminosities previously estimated from HIPPARCOS data, masses from pulsations and C/O abundance ratios from spec- 
troscopy, and effective temperatures from a new homogeneous scale, are used. Predicted and observed mass loss rates agree 
fairly well, as functions of effective temperature. The signature of the mass range M < 4Mq of most carbon-rich AGB stars is 
seen as a flat portion in the diagram of mass loss rate vs. effective temperature. It is flanked by two regions of mass loss rates 
increasing with decreasing effective temperature at nearly constant stellar mass. Four stars with detached shells, i.e. episodic 
strong mass loss, and five cool infrared carbon-rich stars with optically-thick dust shells, have mass loss rates much larger 
than predicted values. The latter (including CW Leo) could be stars of smaller masses (M ^ 1.5 - 2.5 Mq) while M = 4Mq 
is indicated for most of the coolest objects. Among the carbon stars with detached shells, R Scl returned to a predicted level 
(16 times lower) according to recent measurements of the central source. The observed expansion velocities are in agreement 
with the predicted velocities at infinity in a diagram of velocities vs. effective temperature, provided the carbon to oxygen abun- 
dance ratio is 1 < 6c /eo ^ 2, i.e. the range deduced from spectra and model atmospheres of those cool variables. Five stars with 
detached shells display expansion velocities about twice that predicted at their effective temperature. Miras and non-Miras do 
populate the same locus in both diagrams at the present accuracy. The predicted dust-to-gas density ratios are however about 
2.2 times smaller than the values estimated from observations. Recent drift models can contribute to minimize the discrepancy 
since they include more dust. Simple approximate formulae are proposed. 

Key words, stars: AGB and post-AGB - stars: carbon - stars: mass loss 



1. Introduction 

The observations of CO and HCN lines in the millimetric range 
together with infrared fluxes, especially from the IRAS cata- 
logue (1988), provide estimates of parameters of the circum- 
stellar envelopes of carbon-rich giants. They are essentially TP- 
AGB stars, i.e. late stages of stellar evolution of low and inter- 
mediate mass stars, experiencing substantial mass loss. This 
latter phenomenon may terminate their evolution on the AGB 
before rapid nucleosynthesis forces them to leave that branch. 
Mass loss also contributes to the replenishment of the interstel- 
lar medium in carbonaceous material (dust and gas), mostly 
from extreme ("infrared") carbon stars. Models and simplify- 
ing assumptions are required to derive the estimates of the mass 
loss rate M in solar mass per year, the expansion velocity v e 
in kilometer per second, and the dust-to-gas mass density ra- 
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tio. For instance a model for the photodissociation of the cir- 
cumstellar CO by Mamon et al. (1988 1 is usually applied to 
radial brightness distributions with fair success, but clear de- 
viations were noted in a few cases (e.g. Schoier & Olofsson 
2001 1. Predominant species, like hydrogen or helium, are not 
traced out and the adopted abundances influence the estimates 
of the total mass loss rate. The estimated distances are also of 
paramount importance since the rate calculated from observa- 
tions increases as the squared distance. A full discussion is out- 
side the scope of the present paper and we refer the reader to 
the papers cited in Sect. 2. 

It has become clear that the interaction of pulsation and 
dust formation plays a key role in the mass loss phenomenon 
(e.g. Wallerstein & Knapp 1998 for a review). The effects of 
stellar pulsation are simulated in dynamical models by apply- 
ing a piston at the inner boundary. Consistent models of cir- 
cumstellar dust shells around carbon-rich long period variables 
(LPV) include time-dependent hydrodynamics, a detailed treat- 
ment of the processes of formation, growth and evaporation of 
dust grains, a carbon-rich chemistry, and radiative transfer (e. 
g. Fleischer et al. [19921 Fleischer [19941 Winters et al. [1994b 
119951 Hofner et al. 119951 [19 96). The empirical Reimers rela- 
tion for RGB (Reimers ll975' l proved its inadequacy for AGB. 
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Table 1. Mean distances (pc) of carbon-rich giants and dispersions as obtained from data of various authors: 093 stands for 
Olofsson et al. dl993a> . G99 for Groenewegen et al. (JT999), SOI for Schoier & Olofsson J200T1 . G02 for Groenewegen et al. 
(2002bi», and L93 for Loup et al. iTWR . They are compared with astrometric data from Bergeat et al. (2002a 2002b =BKR) for 
(n) stars in common. Ratios to BKR means are then quoted in the last line (See text for a discussion). 
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Alternative relations were derived from observational results 
and/or theoretical models of mass loss (e.g. Volk & Kwok 
[TWgl Vassiliadis & Wood fT^l Blocker [l99ll . The mass loss 
relation of Dominik et al. (1990 1 was the first one based on 
detailed calculations of dust formation and growth (Gail & 
Sedlmayr 1988 1, with rather extreme stellar parameters adapted 
to the very late stages of AGB evolution. Approximate formu- 
lae for mass loss rate, expansion velocity and dust to gas den- 
sity ratio were then calculated by Arndt et al. (119971 . They con- 
sidered six stellar parameters as independent ones, namely the 
temperature, the luminosity, the mass, the abundance ratio of 
carbon to oxygen, the pulsational period and the velocity am- 
plitude of the pulsation. They then derived approximate equa- 
tions on the basis of 48 dynamical models, exploring various 
ranges for the six parameters. They also compared their results 
to those of 22 analogous models from Hofner & Dorfi ( 1997 > 
who proved to be well approximated by their set of equations. 

The predictions of Arndt et al. are compared to the data de- 
duced from observations as published by various authors. We 
use the effective temperatures of Bergeat et al. (2001) as stel- 
lar temperatures. Using a new homogeneous scale constructed 
from model atmospheres, observed angular radii and spec- 
tral energy distributions, they have proved efficient in solving 
various problems. We also adopt the distances and luminosi- 
ties derived by Bergeat et al. {2Q02b i from the HIPPARCOS 
data ( 19971 by taking into account three distinct biases (their 
Sect. 2.3; see also Knapik et al. 1998 and Bergeat et al. 
2002a). Pulsation masses were derived for carbon-rich LPVs 
by Bergeat et al. ( 2002c I. The carbon to oxygen (C/O) ra- 
tios were taken from Lambert et al. (1986 1, Olofsson et al. 
( 1993b I and Abia & Isern ( 1996 1. The periods of pulsation used 
are the photometric periods taken from the General Catalogue 
of Variable Stars (GCVS; Kholopov et al. I1985i or from the 
HIPPARCOS data JT9971 . The values are those of fundamen- 
tal mode pulsators from the study of Bergeat et al. ( 2002c I. 
Considering the complexities of those pulsating atmospheres 
and the uncertainties on transfer in their models, we adopt the 
reference value Au = 2 km s~' of Arndt et al. ( 1997 1 for the ve- 
locity amplitude of the pulsation. We shall see that its influence 
is quite limited. 

We re-scaled the mass loss rates from the literature to the 
distances of Bergeat et al. J2002bt and a f-correction was ap- 
plied to convert them to the scale of Schoier & Olofsson ( 2001 ) 
who properly dealt with radiative transfer (Sect. 2). In addi- 
tion, the expansion velocities and dust-to-gas density ratios 
were compiled for a sample of carbon-rich giants (Sect. 2). The 



available predictions from hydrodynamical models are summa- 
rized in Sect. 3 and equations fitted on a grid from the litera- 
ture are adopted. Values of stellar parameters from Bergeat et 
al. (2001 2002b 2002c) were compiled for stars in common 
with millimetric studies and mean values collected for eight 
photometric groups (Sect. 4) to be replaced in the formulae of 
Arndt et al. ( 1997). Then, the comparison of observed and pre- 
dicted mass loss rates is performed in a diagram against effec- 
tive temperature (Sect. 5). In a similar approach, the observed 
and predicted expansion velocities are compared by plotting vs. 
effective temperature (Sect. 6). The discrepancy noted between 
observed and predicted dust-to-gas density ratios, is estimated 
and discussed (Sect. 7). A few conclusions are given and dis- 
cussed (Sect. 8) and simple approximate formulas are proposed 
in Appendix A. 

2. The data from millimetric observations 

Mass loss rates, expansion velocities and dust-to-gas density 
ratios were extracted from the literature for a large sample 
of carbon-rich giants. The early millimetric data and results 
are summarized in a catalogue by Loup et al. ( 1993). Schoier 
& Olofsson (2001) revised a set of data from Olofsson et al. 
( 1993a 1993b). The mass loss rate, for stars in common, is on 
average 3.1 + 3.6 times larger in the 2001 paper than it was 
in those of 1993. Values in Groenewegen et al. (2002a) for 18 
stars in common with Olofsson et al. are on average 2 times 
larger. A mean 2.9 ratio was found from two stars in common 
with Winters et al. (2003). Admittedly, recent accepted values 
for mass loss rates are typically 2-3 times higher than the pre- 
vious ones from Olofsson et al. ( 1993a 1993b). We also con- 
sidered data for a few stars from Groenewegen et al. (QJ)99 ), 
Schoier & Olofsson (EUUOl and Le Bertre et al. EUOH 
see also Meixner et al. [T998l for RW LMi=CIT6 and Skinner et 
al. ll998l for CW Leo=IRC+10216). 

Usually, the distances adopted by the authors substantially 
differ and this is a fundamental problem since the mass loss rate 
(for both gas and dust) increases as the squared distance. Large 
discrepancies may exist between distances estimated from var- 
ious methods: mean absolute magnitude adopted, say in the 
K-band, mean bolometric magnitude with bolometric correc- 
tion or mean luminosity, absolute magnitude vs. period for 
long period variables, kinematical methods ... The astromet- 
ric parallaxes of those bright distant giants are also affected 
by substantial errors. First of all, we calculated the mean dis- 
tances and dispersions for the various samples of the above- 
mentioned references and compared them to their counterparts 
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Table 2. The data for 119 carbon-rich giants (Sect. 2) with variable star names (Columns 1 and 10, except for CGCS1006 in 
Stephenson 1989 and Alksnis et al. 2001 1 shown with boldface characters for 19 Miras. Data from SEDs and spectra are given: the 
photometric groups in Columns 2 and 1 1 with an additional J for 13 C-rich stars, the effective temperatures (Kelvin) in Columns 
3 and 12, the carbon to oxygen abundance ratio (C/O) of the atmosphere in Columns 4 and 13, the absolute bolometric magnitude 
M[, in Columns 5 and 14, that can be translated into luminosity in solar units from Eq. (3), and the distance D BKR in parsec in 
Columns 6 and 15. The data from millimetric observations are displayed in Columns 7 and 16 (mass loss rate as M' = 10 8 M 
where M in solar mass per year), Columns 8 and 1 7 (expansion or terminal velocity v e in km s~ 1 , and Column 9 and 1 8 (r ' = 1 3 r 
where r = Md/M g is the dust-to-gas ratio of mass loss rates). Additional notes deal with ranges: a HC5-CV4 & 3525-2775 K, b 
HC5-CV6 & 3520-2645 K & M b = -3.72; -4.36, c 2290-2245 K&M b = -4.86; -4.53, d 2920-2840 K & M b = -5.22; -4.80, 
e M b = -5.81; -5.48, f 3345-2825 K & M b = -6.15; -5.67, e 2405-2525 K & M b = -5.89; -6.10, h IRC+10216 1915-2105 K, 
'CIT6 2425-2465 K, J 2100-1945 K & M b = -5.23; -4.60, k 2865-2680 K & M b = -5.51; -5.00, 1 CV5-CV6 & 2650-2530 K & 
M b = -4.94; -5.06, ra 1885-1875 K & M b = -6.28; -5.93, n 2945-2865 K & M b = -4.84; -5.61, 1975-1875 K, p 2240-2095 K 
&M b = -5.42; -5.22, q CV1-CV5 & 3245-2605 K & M b = -3.01; -3.69, 1 2735-2655 K & M b = -4.78; -4.58, s CIT 5, 1 IRAS 
05352+2247, u AFGL 971, v AFGL 2067=IRC-10396, w AFGL 31 16=IRC+40540, x AFGL 1235, y AFGL 1961 (erroneously 
identified as V5220ph). Not shown in Fig.[2: according to Olofsson et al. (2000 1, TT Cyg (which is surrounded by a thin shell) 
is presently losing mass at a modest rate of 3 10 _8 Mq yr~' which is perhaps to be multiplied by the 3.9 correction factor of Eq. 
(2); according to Izumiura et al. ( 1996 1, Y CVn** exhibited a 7 - 20 10~ 6 Mq yr~' rate at the formation of its thin shell. 
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from Bergeat et al. J2002bl BKR) whose values are astromet- 
ric in essence and deduced from HIPPARCOS data (Knapik et 
al. 1998 and Bergeat et al. 2002ai, with no photometric criteria 
used. The results are shown in Table 1 . There is apparently an 
excellent agreement for the 104 stars of Olofsson et al. ( 1993ai 
in common, with 707-708 pc as mean values. A detailed anal- 
ysis (Bergeat 12004ft however proved this is fortuitous since a 
Malmquist-type bias occurred in the Olofsson et al. distances 
in this sample of magnitudes K < 2, adopting Mk - -8.1 (ac- 
tually there is a 3-magnitude range whose central value is close 
to -7.50). Conversely, no such bias appeared when compar- 
ing the distances deduced from the Mic-magnitudes derived 
by Mennessier et al. (2001 1 through a maximum likelihood 
method for 124 carbon giants in common with BKR. There is 
also a good agreement with BKR in Table 1 for a small sample 
of 11 stars from Groenewegen et al. ( 1999 1. They are located 
farther away at 1141-1161 pc on average, that is 1.63 times 
more distant. We also considered a sample of 54 stars in com- 
mon to Olofsson et al. ( 1993a) and Schoier & Olofsson ( 2001 1. 
There is a 1.33 ratio (i.e. 1.76 on mass loss rates), in the di- 
rection opposite to that of the 3.1 mean ratio of mass loss rates 
mentioned above. It is to be compared to 0.92 when referred to 
BKR. The Groenewegen et al. (I2002b> mean distance is 2.05 
times larger than the BKR-value for 25 stars in common, lead- 
ing to their mass loss rates being systematically larger. The fac- 
tor is only 1.35 between Loup et al. <!1993t and BKR. 

Finally, for the sake of consistency, we transformed the 
rates M given at distances D in the mentioned papers, accord- 
ing to the relation 

M BKR = M (Dbkr/D) 2 (1) 

where Dbkr stands for the distances of Bergeat et al. ( 2002b I. 
The gas mass loss rates were derived from CO observations 
by a formula established by Knapp & Morris ( 1985 1 for unre- 
solved optically thick envelopes. From a more detailed radia- 
tive transfer analysis, Schoier & Olofsson (1200 1> found that 
the rates resulting from the formula are systematically under- 
estimated when compared to the new values (their Fig. 8). As 
expected, the discrepancy increases with decreasing rate i.e. de- 
creasing optical depth. From the comparison of Mbkr for 53 
stars in common we obtained a mean correction factor (f) of 
2.5 + 1.7 to scale the rates from Olofsson et al. ( 1993ai to the 
Schoier & Olofsson (2001 1 level. We adopted mean values in 
four ranges 

0.94 + 0.34 1 if M BK R > 4 10~ 6 
2.2 + 1.0 if 4 10~ 7 < M BK R < 4 10~ 6 
2.5 + 1.4 if 7 10~ 8 < M BK R < 4 10~ 7 

3.9 + 3.1 if M BK R < 7 10~ 8 (2) 

where quoted rates refer to values from Olofsson et al. 
(1993a|. This is in fair agreement with Fig. 8 in Schoier 
& Olofsson (2001). About 70% of the data we used were 
taken from those two references. From 17 stars in common, 
we found f = 1.9 + 1.4 for Groenewegen et al. ( I2002bl 8%) 
while f= 3.3 + 0.3 was deduced from the seven values of 
Groenewegen et al. (1999 2.7%). Here too, there is some in- 
dication of f tending to unity for objects with very large rates 



but we were unable to split the samples. For Loup et al. ( 1993 
7 stars and 9%), we obtained f = 1.2 ± 0.5 while f = 3.4 + 2.1 
and f = 3.1 + 1.7 were deduced for Le Bertre et al. (2001 5 
stars) and ( 2003 4 stars) respectively. These later five correc- 
tion factors and those in Eq. (2) were applied to the MBKR-data 
derived from Eq. (1). We rejected the data available in Winters 
et al. ( 2003 1 since the correction factor curiously amounted to 
0.4 from only two stars in common with Schoier & Olofsson. 

Both the obtained mean gas mass loss rates and the 
DBKR-distances are given in Table 2 together with photometric 
groups and effective temperatures from Bergeat et al. (|2001 1. 
The luminosities in solar units can be derived from the abso- 
lute bolometric magnitudes of Bergeat et al. ( I2002b> making 
use of 

L/L = i -°- 4 ^- 4 - 73 ) . (3) 

Additional star identifiers and detailed bibliographic sources of 
the mass loss data for each star are provided in Table 3, only 
available in electronic form. The C/O ratios quoted in Table 2 
were taken from Lambert et al. ([1986 1, Olofsson et al. ([1993b I 
and Abia & Isern ( 1996). Also given are the expansion veloc- 
ities from the same five references mentioned above for mass 
loss rates. Both the dust and gas mass loss rates are proportional 
to the adopted squared distance and no corresponding correc- 
tions need be applied to the dust-to-gas ratios. The dust mass 
loss rate is derived from the IRAS (60^m) making use of a for- 
mula given by Sopka et al. ( 1985] see e.g. Eq. (8) in Olofsson 
et al. ( 1993ai. The dust rate is proportional to L~ . We thus 
scaled the values from Olofsson et al. ( 1993b), Groenewegen 
et al. ( 1999) and Groenewegen et al. ( 2002b I to the luminosi- 
ties quoted in Table 2, according to the individual luminosities 
adopted by these authors. The obtained values were then di- 
vided by the correction factors of Eq. (2) or by 3.3 or 1.9, i.e. 
the f-factors for the gas mass loss rates from those three refer- 
ences. As discussed in Sect. 7, the latter correction strongly 
shifts observations toward model predictions. Two estimates 
are provided by Olofsson et al. ( 1993a), namely a lower one as- 
suming no drift velocity for the dust, and an upper one adopting 
an upper limit for drift velocities (up to 20 km/s, i.e. eventu- 
ally higher than the gas expansion velocity) obtained by equat- 
ing the radiation force on the grains and the drag force due 
to gas-grain collisions in the limit of supersonic motion (e.g. 
Kwok fT975l . The estimates from Groenewegen et al. make use 
of moderate drift velocities, usually 2-3 km/s. Their mean val- 
ues, (1.2 + 1.3) 10 3 for Groenewegen et al. Il999l n=8) and 
(1.45 ± 1.1) 10~ 3 for Groenewegen et al. ( I2002bl n=24) are 
consistent with that of the "lower" values ((1.3 + 1.0) 10~ 3 , 
n=55) of Olofsson et al. ( 1993a), while the "upper" values yield 
a significantly larger mean ((3.2 + 2.5) 10~ 3 , n=55). Average 
values of the available data are quoted in Table 2, ignoring the 
"upper" ones in Olofsson et al. fT993at. They thus correspond 
to the low drift velocities adopted here. 

The sample of Table 2 has 119 carbon-rich giants with at 
least one of the three quantities from millimetric data given. 
Many "optical" carbon stars of low to moderate mass loss rates 
are included while the "infrared" carbon stars with higher rates 
and often optically thick dust shells are fewer. This is due to 
our requirement that stellar data be available from Bergeat et 
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al. ( EUTTTI I2002bl I2002c> that, in turn, necessitated sufficient 
photometric data in the visible range. This is very often miss- 
ing for "infrared" carbon stars. This is rather unfortunate since 
330 "infrared" carbon stars were studied by Groenewegen et al. 
j2002all2002bl 

3. The predictions from hydrodynamical models 

As mentioned in Sect. 1, both pulsation and dust formation 
are involved in the mass loss phenomenon on the asymptotic 
branch. The mass loss relation of Dominik et al. ( 1990) is based 
on calculated models and was the first one available. It de- 
scribes the stationary outflows of carbon-rich low mass stars. 
The mass loss rate increases with increasing luminosity and 
decreasing effective temperature. The strong dependence upon 
the effective temperature is due to the fact that dust always con- 
denses at about the same temperature. For stars with higher ef- 
fective temperatures, this condensation temperature is reached 
farther out in the shell where dust formation is less effective. 
The mass loss rate increases as the stellar mass decreases. 
Other things being equal, increased gravity due to larger masses 
inhibits the wind. The dependence of the mass loss rate on the 
carbon to oxygen ratio is much less marked, especially for val- 
ues larger than 1.5. Dominik et al. proposed an analytical fit, 
their formula (7), which is presumably accurate for extreme 
mass loss rates, i.e. for M > 10~ 6 MQyr~'. The terminal veloc- 
ity ranges from 5 to 40kms~ I . It increases with decreasing ef- 
fective temperature and with increasing luminosity or increas- 
ing C/O ratio. The authors proposed a linear fit as their Eq. (8). 
The dust-to-gas density ratio is closely related to the final de- 
gree of condensation, as shown by formula (9) of Dominik et 
al., and it ranges from 6 10" 4 to 5 10~ 3 . 

Making use of the models of Fleischer et al. J1992I and 
Fleischer (IT9941 . Arndt et al. ( fT997t derived a grid of 48 dy- 
namical models usually with less extreme mass loss than in 
the cases explored by Dominik et al. (I1990> . They are real- 
istic models of circumstellar dust shells around carbon-rich 
long-period variables. They include the calculation of time- 
dependent hydrodynamics as well as a detailed treatment of 
dust formation, growth and evaporation. The maximum likeli- 
hood method was applied to the grid with six parameters they 
considered as independent, namely the stellar temperature, the 
stellar luminosity, the stellar mass, the abundance ratio of car- 
bon to oxygen, the pulsational period and the velocity ampli- 
tude of the pulsation. They obtained their equations (1), (3) and 

(5) for the mass loss rate, the dust-to-gas ratio and the expan- 
sion (outflow) velocity, respectively. The trends are essentially 
the same as those obtained by Dominik et al. and summarized 
above. Keeping constant the parameters of low influence, they 
also wrote the simplified relations they numbered (2), (4) and 

(6) for the same three "mass loss" quantities. 

Applying a piston and a variable luminosity at the inner 
boundary in the models of Hofner et al. ([1995 1996 1, Hofner 
& Dorfi ( 1997 1 constructed a set of 22 carbon-rich LPV mod- 
els. The ranges of the parameters were smaller but Arndt et al. 
( 1997 1 conclude a good agreement with their work for values 
in common. 



We recall here the equations (1), (3) and (5) from Arndt 
et al. ( 1997 1 whose predictions we intend to compare with the 
observed data as collected in Sect. 2. For the mass loss rate in 
solar mass per year they derived 

logM = -4.95 - 9.45 log (T/2600) + 1.65 log (l/1O 4 L ) 
-2.86 log (M/M©) + 0.47 log (e c /l .8e ) 
-0.146 log (P/650) + 0.449 log (Au/2.0) (4) 

where the parameters of major influence are the effective tem- 
perature T in Kelvin, the luminosity L and the mass M. The mi- 
nor parameters are the carbon to oxygen abundance ratio, the 
period of the LPV in days, and the piston velocity amplitude in 
km s _1 . For the dust-to-gas density ratio, they obtained 

log (p d /p g ) = 2.74- 1.62 logT- 0.248 log (L/L©) 
-0.658 log (M/M ) +2.45 log(e c /e ) 

+0.230 log P -4.93 10~ 3 log (Au) (5) 

with the same units. In addition to T, L and M, the carbon 
to oxygen abundance ratio is here a major parameter. Again, 
the period and piston velocity amplitude have little influence. 
Finally, they give for the expansion (or outflow) velocity 

log Voo = 0.730 - 0.148 logT + 0.107 log (L/L Q ) 
-6.40 10~ 2 log (M/M ) + 1 .74 log (e c /eo) 

+9.56 10~ 2 log P - 7.07 10~ 4 log (Au) (6) 

where the effective temperature, the luminosity and the carbon 
to oxygen abundance ratio prove to be the main parameters. 
The stellar mass, the period of pulsation and the piston velocity 
amplitude have little influence here. 

4. The realm of carbon-rich AGB stars 

The adopted values of the six parameters in the models men- 
tioned in Sect. 3 were varied more or less freely by the authors. 
The resulting Eqs. (4), (5) and (6) are given for ranges of pa- 
rameter values that are those of the used (satisfactory) mod- 
els. Each parameter is supposed to vary independently which 
is actually not the case for existing carbon-rich giants on the 
AGB. We briefly summarize here the results from Bergeat et 
al. (p001 2002b 2002c I and references therein, on this spe- 
cific topic. The bright carbon-rich giants are located on the 
asymptotic branch, with various initial masses and chemical 
compositions. Consequently, those Galactic giants populate a 
strip in the HR diagram with increasing luminosity for de- 
creasing effective temperature (See Figs. 8 and 9 of Bergeat 
et al. 2002b I. Mean values of both quantities were given in 
their Table 3 for the HC-CV groups. This classification was 
shown to be strongly correlated with effective temperature de- 
creasing along the sequence HC0 to HC5 and then CV1 to CV7 
(Bergeat et al. 2001 1. Influenced by true dispersion and errors 
on distance estimates, the populated strip is rather large but the 
trend in mean values is obvious. A study of pulsation modes 
resulted in the evaluation of pulsation masses in good agree- 
ment with evolution tracks in the HR diagram (Bergeat et al. 
2002c I for long period variables. Mean pulsation masses and 
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periods increase with decreasing effective temperature along 
the sequence of photometric groups, as shown in their Table 
1. In a mass-luminosity diagram (their Fig. 4), the mean lumi- 
nosity increases with increasing mean mass. For the cool car- 
bon variables (CV-stars), the carbon to oxygen abundance ratio 
was plotted against effective temperature in Fig. 7 of Bergeat 
et al. (2002c). The mean values and dispersions are increasing 
with decreasing effective temperatures along the CV1 to CV6 
sequence. They however decrease below 2500 K at the CV6- 
CV7 junction, a result possibly due to large dust condensation 
(SiC and carbonaceous grains) in the atmosphere and resulting 
carbon depletion in the gas phase. The only parameter of for- 
mulae (4) to (6) in Sect. 3 for which we have no indication, 
is the piston amplitude velocity. It is however supposed to be 
higher in Miras and semi-regulars with large amplitudes than 
it is in semi-regulars with low amplitudes and irregular vari- 
ables. Its influence is however small on the three quantities cal- 
culated from Eqs. (4) to (6). We thus adopted the Au = 2 km s _1 
value as favored by Arndt et al. ( 1997 1. We have summarized 
the values in Columns 7 to 9 of Table 4 as calculated from the 
mean data of Columns 2 to 6, following the same sequence of 
photometric groups. In those calculations, we eventually went 
beyond the validity ranges of formulae (4) to (6) stricto sensu, 
as defined by the parameter values that Arndt et al. introduced 
in their 48 models. Discussions are postponed to the following 
sections. 

5. Mass loss rate vs. effective temperature 

The mass loss rates M as quoted in Table 2 (Columns 7 and 16) 
are plotted on a logarithmic scale in Fig.[Qas a function of the 
effective temperatures from Columns 3 and 12. The main trend 
of increasing mass loss rate with decreasing effective temper- 
ature is observed as expected from theory (Sect. 3). Typical 
error bars on both coordinates are also shown. The Miras are 
displayed as crosses while the filled circles stand for semi- 
regulars. Among the latter, four stars with detached shell(s) 
are shown with x-symbols in the 2600-2800 K range. Denoted 
by asterisks in Table 2, they are U Ant, S Set, R Scl and U 
Cam. The observation of detached shell(s) is indicative of re- 
cent enhanced mass loss. As a consequence, those four stars 
do exhibit larger mass loss rates than the bulk of our sam- 
ple, in the same temperature range. Recent observations of a 
HCN line emission in the central source of R Scl (2625 K) that 
was resolved were secured with a FWHM of 1" by Wong et 
al. (|2004i. The authors thus derive for the present mass loss 
rate 2 lO~ 7 M0yr _1 for an adopted distance of 360 pc, a rate 
we convert into 3.5 lO^Mgyr -1 for our 475 pc distance esti- 
mate (see Column 15 of Table 2). This latter value, displayed 
as a diamond-symbol in Fig. ^ is located in the main locus. 
This is about 16 times less than the earlier 5.5 lO^M^yr" 1 
value quoted for R Scl in Column 16 of Table 2, a rather drastic 
change. This is convincing evidence for highly episodic mass 
loss, possibly caused by a thermal pulse. The detached shell of 
R Scl has recently been investigated for MgS, mining the mass- 
loss history of this object (Hony & Bouwman 1998 see also 
Olofsson et al. 2000 for a detailed study of circumstellar CO 
in TT Cyg). We conclude that mass loss rates larger than about 



Table 4. The mean stellar data for carbon-rich giants (Columns 
2-6) from Bergeat et al. d200T1 12002bl I2002c> and the corre- 
sponding mass loss data (Columns 7-9) calculated from for- 
mulae (4), (5) and (6) taken from Arndt et al. J1997i . For 
photometric groups of Column 1, the mean values of effec- 
tive temperatures in Kelvin (Column 2), the mean luminosi- 
ties 1 =< L/Lq > (Column 3) and masses m =< M/Mq > 
(Column 4) in solar units, the mean carbon to oxygen abun- 
dance ratio C/O (Column 5) and the mean pulsation period P in 
days (Column 6). The group CV5 "und" corresponds to a small 
sample of underluminous CV5-stars which are not included. 
Adopting the piston velocity amplitude Au = 2kms~' favored 
by Arndt et al., the corresponding mass loss rates M' = 10 7 M 
(Column 7) where M is in solar mass per year, the dust-to-gas 
abundance ratios p<j/ g = 10 4 pd/Pg (Column 8) and the expan- 
sion velocities v e in km s (Column 9). 
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10~ 6 M©yr~' in the 2500-2800 K range are representative of 
temporarily-enhanced mass loss leading to detached shell(s) 
hundreds or a few thousands of years later (see e.g. Lindqvist 
et al ll996l for U Cam, Izumiura et al. 1997 for U Ant, Izumiura 
et al. lT996l and Le Bertre & Gerard l2"Urj4l for Y CVn). 

At lower temperatures, we used symbols surrounded by 
squares for five stars with optically-thick dust shells as shown 
by large infrared excesses. They are CW Leo (IRC+10216), V 
346 Pup and LP And (three Miras), and RW LMi (CIT 6) and 
FX Ser (two semi-regulars). They exhibit mass loss rates much 
larger than those from the carbon-rich giants with optically- 
thin dust shells in the same 2500-1800 K range, that is larger 
than about 6 lO^MQyr -1 for effective temperatures less than 
2500 K. A few intermediate objects are also noticed (V Hya, 
BN Mon). The circumstellar envelope characteristics of 9 car- 
bon stars were probed by Schoier et al. ( 2002 ), combining radio 
measurements and infrared observations from ISO. Their radia- 
tive transfer analysis allowed us to put constraints on the mass 
loss history of those objects, specially CW Leo=IRC+10216 
and RW LMi=CIT6 (studied here) and IRAS 15194-5115=11 
Lup=CGCS 3592 (not included). They found that those stars 
are not likely to have experienced any drastic long-term mass 
loss rate modulations, at least less than a factor of about 5, over 
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Fig. 1. Logarithm of the mass loss rate in solar mass per year as a function of effective temperature in Kelvin, as taken from 
Table 2 (observations). Typical error bars are displayed at left. Semi-regulars and Miras are shown as filled circles and crosses 
respectively with 4 non-Miras with detached shells (x and diamond symbol for present R Scl), and 5 objects with optically- 
thick dust shells (squared symbols). For a few variables, there are two symbols for temperatures obtained at distinct phases. The 
predictions from models (Table 4) are plotted as filled squares connected by full lines. Additional predictions for various mean 
masses are shown as filled triangles connected with dashed lines (see text). 



the past thousands of years. This is a situation in contrast with 
the case of the detached shells mentioned above. 

If we remove the nine symbols of those high mass loss ob- 
jects in Fig. [2 the remaining filled circles (semi-regulars) and 
crosses (Miras) populate a locus with a flat portion in the 2900- 
2400 K range and for n=66 stars, 

logM (-0.8 + 1.7) log T eff - (3.6 + 5.8) (7) 
with an average value of 

< M >^ (5.8 + 4.4) 10~ 7 M e yi-- 1 (8) 

for this wide 10~ 7 - 2 10~ 6 MQyr~' nearly horizontal strip. It 
is flanked by two ranges of increasing mass loss rate with de- 
creasing effective temperature. For T e ff > 2900 K we obtain for 
32 stars 

M ^ 4.51 10 2() T;™ 3 Mqyt- 1 (9) 
and for T eff < 2400 K with 17 stars 

M^^IO 24 ^ 20 M yr-' . (10) 

Those three relations which are not satisfied by the 9 objects 
with strong mass loss and possibly a few others, are not shown 
in Fig. [2 Instead the predictions of Sect. 4 as quoted in Table 4 
are plotted as filled square-symbols connected by a full broken 
line. A very good agreement is observed in the 3000-2400 K 



range of effective temperatures and the observed flat portion is 
confirmed by the predictions. We shall identify its cause later 
in this section. Below 2400 K, the observed locus is well repro- 
duced by predictions provided M = 4.2 M0 is adopted instead 
of the uncertain M = 8.2 Mq value estimated by Bergeat et al. 
(2002c 1 from a few extreme objects. This high value was pos- 
sibly induced by overestimating the "photospheric" radius of 
very cool giants with dust shells. The right portion of Fig. ^ 
is an additional argument in favor of masses not exceeding 
about 4 Mo for a large majority of luminous carbon-rich gi- 
ants. For effective temperatures higher than 3000 K, there is an 
increasing disagreement between observations and predictions 
up to 3400 K. Hot variables exhibit an irregular pulsational be- 
haviour with small amplitudes, resulting in dubious pulsation 
mode and period. Uncertainties increase accordingly. It should 
be noted that, according to Arndt et al. (fT9"97). one has to re- 
frain from applying Eqs. (4) to (6) beyond 3000 K since none 
of their 48 models lies there. In addition, none of their mod- 
els is fainter than 5 10 3 L© contrary to the observed averaged 
mean luminosities of HC5, CV1 and CV2-stars in Table 4. We 
however attempted to apply the < M >= 1.6 Mq mean value 
of the CV3-group in Table 4 to the other earlier groups CV1 
and CV2 with other data otherwise unchanged. Together with 
an additional calculation from Eq. (4) assuming for CV1 the 
< M >= 1 .0 M© value quoted for CV2 in Table 3, they are plot- 
ted as filled triangles in Fig.^ The dashed lines connecting the 
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Fig. 2. Outflow (expansion) velocities in km/s as a function of effective temperatures in Kelvin, as taken from Table 2 (observa- 
tions) and Table 4 (predictions from models). Typical error bars are also displayed at left. Same symbols used as in Fig.^ with 
5 non-Miras with detached shells (x) . The predictions from theory are plotted as filled squares connected by full lines. Extreme 
lines corresponding to C/0=1 and C/0=2 respectively are also shown (see text). 



CV3-square to the filled triangles for CV1 and CV2 and the 
initial segmented full line reasonably well encompass the ob- 
servations. This result suggests that the slope for T e fr > 2900 K 
in Fig. [2 is indicative of a nearly uniform value of about one 
solar mass. A large correction factor (3.9 or 2.5) from Eq. 
(2) was applied to part of the contributing data. Without that 
correction, the shift of observations below the continuous line 
of predictions would have been substantially larger. Similarly 
the slope below 2400 K seems to correspond to a nearly uni- 
form value of about 4Mq. The latter is predicted as the up- 
per limit from evolutionary calculations presumably due to hot 
bottom-burning most effective in AGB stars of larger masses 
(Bergeat et al. 2002b and references therein). We also calcu- 
lated the mass loss rate of the CV6 and CV7-groups, assuming 
the < M >= 2.6 M0 value quoted for CV5 in Table 3, keeping 
the other parameters unchanged. They are shown as filled tri- 
angles connected to the CV5-square by a broken dashed line. A 
possible explanation of higher mass loss rate might be smaller 
mass due to previous strong mass loss. Precisely, Schoi'er et 
al. (|2002 1 deduced that mass loss rates remained high over the 
past thousands of years in objects like CW Leo=IRC+10216 
and RW LMi=CIT6. From the comparison of theoretical nu- 
cleosynthesis models and measurements of abundances in the 
circumstellar envelope, Kahane et al. (2000) favor a mass not 
larger than 2M© for CW Leo, well below the 4.2M0-limit 
considered for the stars of Eq. (10). Winters et al. ( 1994ai pre- 
viously estimated M < 2M with C/O 1.4 for CW Leo. In 
other words, the apparent (lasting) "superwinds" could well be 



standard mass loss in objects of diminished masses. It should 
be noted that Arndt et al. ( 1997 1 did not calculate models cooler 
than 2300 K. Additional data and models are clearly necessary 
before a firm conclusion can be reached. 

We shall see in Sect. 6 that the four stars with detached 
shells mentioned above also have expansion velocities larger 
than average. Actually, it happens that velocities are little in- 
fluenced by adopted masses (Sects. 3 and 6). In addition, R Scl 
has now returned to low mass loss rates (Wong et al. 2004l. For 
those stars with detached shells and intermediate (2400-2900 
K) effective temperature, the above explanation of enhanced 
mass loss in terms of diminished masses is not acceptable. By 
continuity, we are lead to the conclusion that the 3000-2400 
K flat portion corresponds to the approximate 1-4 Mq mass 
range. Discussing their formula (1), i.e. our Eq. (4), Arndt et 
al. (I1997> concluded that the abundance ratio, pulsation period 
and piston amplitude velocity have only small influences on 
the mass loss rate. Accordingly, we find no clear difference be- 
tween the loci of Miras and semi-regulars. It remains to con- 
sider the effective temperature, the luminosity and the mass of 
the star. Since, on average, luminosity increases with decreas- 
ing effective temperature, both parameters contribute to an in- 
crease of the mass loss rate. Following Eq. (4), that increase is 
inhibited in the 3000-2400 K range by the increase of stellar 
mass in the 1-4 Mq domain. The 3000-2400 Kflat portion in 
Fig-U\is a further argument in favor of the 1 -4Mq range for 
bright carbon-rich giants (less luminous objects can populate 
the 0.5 - 1 Mq domain), in addition to the use of evolutionary 
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tracks in the HR diagram and to the calculations of present pul- 
sation masses (Bergeat et al. 2002b 2002c). We conclude that 
the models and relations from Arndt et al. ( 1997) satisfactorily 
reproduce the mass loss rate deduced from millimetric observa- 
tions. Larger mean mass thus implies, other things being equal, 
a lower mass loss rate. 



6. Outflow velocity vs. effective temperature 

The outflow or expansion velocities v e , i.e. estimates of veloc- 
ity at infinity in km/s as quoted in Table 2 (Columns 8 and 17) 
are plotted in Fig.[2]as a function of the effective temperatures 
from Columns 3 and 12. The main trend of increasing expan- 
sion velocity with decreasing effective temperature is observed 
as expected from theory (Sect. 3), but with a large scatter here. 
Typical error bars on both coordinates are also shown and are 
much smaller. The same symbols as in Fig. ^ are used. The 
X-symbols point to five stars with detached shells (U Ant, S 
Set, R Scl, U Cam and TT Cyg). The latter also show large ve- 
locities (14-21 km/s) for T eff 2600 - 2800 K about a factor 2 
larger than predicted values. We already found in Sect. 5 that 
they display higher mass loss rates than predicted. The issue is 
dubious for the five cool extreme objects of Sect. 5. They are 
scattered with CW Leo=IRC+10216 on the predicted relation 
and FX Ser located well above it. 

The predictions of v M from Eq. (6) as quoted in Table 4, are 
shown in Fig.fJ]as filled square-symbols connected by a broken 
full line. The latter fits reasonably well the 86 % of the stars in 
the sample that are located in the lower part of the diagram. The 
models seem to correctly predict the low expansion velocities 
which roughly correspond to the low mass loss rates of Sect. 5, 
themselves well-predicted by the same models. The decreasing 
effective temperature and increasing luminosity along the se- 
quence of Table 3 both contribute to the increase of the outflow 
velocity but the effect is small when compared to the influence 
of the ec/eo abundance ratio. This point is illustrated by the 
two dashed lines in Fig.[5]as obtained keeping the ratio to 1 or 
2 instead of the mean values of Table 4. Within errors, the ob- 
served values lie between those two lines. This is precisely the 
range of the abundance ratio as derived from spectra and model 
atmospheres of the cool carbon-rich giants (Lambert et al. 1986 
and references cited in Sect. 2; see also Fig. 7 of Bergeat et al. 

We also note that above 2900 K the observations in Fig. [2] 
tend to concentrate towards the C/0=1 line, below 10-15 km/s. 
This can be attributed to the fact that stars in the hotter domain 
of Table 4 do have mean abundance ratios close to unity with 
a small range apart. Below 2900 K (CV3 and later groups), 
the C/O mean values and ranges increase (Table 4). On av- 
erage, larger velocities are reached there. The influence of 
the C/O ratio is responsible for a triangular void in Fig. |2] at 
T eff > 2800 K that is inserted between the C/0=1 and C/0=2 
lines. Predictions are consistent with millimetric observations 
and visible and infrared spectroscopy with I < ec I eo < 2, and 
with the mean values and dispersions as functions of effective 
temperature from Bergeat et al. (2001), as displayed in Fig. 7 
of Bergeat et al. [2002c). 
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7. Dust-to-gas ratio vs. effective temperature 

The diagram of the dust-to-gas ratio against effective tempera- 
ture is not shown here. There is actually no marked slope when 
Eq. (5) and the corresponding predictions of Table 4 are used. 
In the 2300-3000 K validity domain of models from Arndt et 
al. dl997> . the mean value amounts to 

< Pd/g >=< p d /p g (5.9 ± 0.7) 10- 4 (11) 

As described in Sect. 2, the mean observed values quoted in 
Table 2 are typically "lower" ones derived assuming low drift 
velocities. For 72 documented stars, we obtained 

M d /M g (1.3 ± 1.1) 10~ 3 (12) 

a value about 2.2 times larger than the predicted one. There is 
practically no slope in the diagram of "observed" dust-to-gas 
ratios of rates against effective temperatures from Table 2, at 
least within the large error bars. It must be emphasized here 
that the data used was corrected for the f-factor of Sect. 2 and 
for a smaller factor for adopted luminosity 

(M d /M g ) BKR = (Md/M g ) L . t (L BKR /L Ll t)- a5 (13) 

Without that f-correction, the above discrepancy ratio would 
have jumped from 2.2 to about 6. The mean "upper" value 
of (3.2 ± 2.5) 10~ 3 derived in Sect. 2 is about 5.4 times larger 
than the Eq. (11) predicted value. Apparently, adopting larger 
drift velocities could make the discrepancy between observa- 
tions and predictions even worse, but predicted values would 
also be increased. Very recent models with the influence of the 
grain drift velocity properly treated (Sandin & H6fner|2004i 
can contribute to shift predictions closer to observations. The 
derivation of the ratio from the observations makes use of data 
from both infrared and millimetric ranges, with simpliflying 
assumptions (like a unique temperature for instance). It is not 
clear whether the same layers are involved in the extended at- 
mosphere and shell (from inner dust condensation zone to outer 
regions with terminal velocity reached). 

Following Eq. (5) taken from Arndt et al. ( 1997 1, it appears 
that in addition to decreasing effective temperature and increas- 
ing luminosity whose influences counteract here, the stellar 
mass and the ec/eo abundance ratio do have an appreciable in- 
fluence. The predicted dust-to-gas density ratio increases with 
both decreasing stellar mass and increasing abundance ratio. 
Adopting for instance M = 1 M0 and ec/eo > 2 can rise the 
Eq. (11) value higher than that in Eq. (12), reaching the level of 
the highest "lower" values observed. It seems outside the range 
of acceptable parameter values taking into account the analyses 
of Sects. 5 and 6 from Table 4 data. However, keeping the mass 
range of Table 4 that fits nicely the gas mass loss rate data (Sect. 
5), an increase of about 20% of the C/O ratios from Lambert et 
al. i 19868 ) can help, which remains compatible with the veloc- 
ity comparison in Fig. [2] Their abundance ratios for 30 carbon 
stars lead to < ( 12 C + 13 c) /H >^ (5.5 + 1.6) 10~ 4 in numbers 
of atoms and (2.7 + 1 .0) 10~ 3 in mass for a gas principally of H 
and He, if their C/O ratios are kept. The estimate is then of 22% 
of carbon locked up in grains if prediction of Eq. (1 1) is trusted 
against 48% for the mean observed ratio of Eq. (12). If con- 
firmed, such figures would mean that carbon molecules could 
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be appreciably depressed. Very high drift velocities could lead 
to unreasonably high estimates. 

8. Discussion and conclusions 

We have compared the results of millimetric (and infrared) ob- 
servations, namely gas mass loss rate, expansion (outflow) ve- 
locity and dust-to-gas density ratio, to the predictions from hy- 
drodynamical time-dependent models of dust-driven winds for 
carbon-rich chemistry. Corrections were applied to the pub- 
lished "observed" values of the first (distance and radiative 
transfer) and the third quantity (luminosity and radiative trans- 
fer) as described in Sect. 2. The prediction of models are ap- 
plied to long period variables whose pulsation is simulated by a 
piston of given amplitude velocity and period acting at the basis 
of the atmosphere. We made use of the approximate formulae 
of Arndt et al. (1997) based on 48 models exploring ranges 
in the six parameters they include. Among them, the pulsa- 
tion period and the piston amplitude velocity have only little 
influence (Sect. 3). We found no marked difference in the ob- 
served data between Miras and non-Miras (here semi-regulars 
with available periods). Four parameters remain whose ranges 
may influence the results, namely effective temperature taken 
as the stellar temperature used in the models, luminosity, mass 
and carbon to oxygen abundance ratio. The observational data 
was summarized in Table 2 for 119 carbon-rich giants formerly 
studied by Bergeat et al. ( 200T1 I2002bl and l2002c> . for effective 
temperatures, luminosities, stellar masses and C/O abundance 
ratio. The formulae (4), (5) and (6) as taken from Arndt et al. 
lead for the mean values of Bergeat et al. to the estimated means 
of Table 4 for the samples corresponding to the photometric 
groups HC5 and C V 1 to C V7 with a range of effective temper- 
atures from 3470 K down to 1950 K. This latter parameter is 
considered as the most accurate one at present and it is used 
for reference. We also hold fixed to reference values the stellar 
mass and the C/O abundance ratio keeping unchanged the run 
of effective temperatures and luminosities of Table 4, namely 
the TP-AGB of Galactic carbon-rich giants in the HR diagram, 
as obtained by Bergeat et al. ( 2002c I. 

As described in Sect. 5, we obtained a good agreement be- 
tween predictions and observations in the diagram of mass loss 
rate vs. effective temperature (Fig. For about 90% of our 
sample, the mechanism involving gas lifted by pulsation and 
radiation pressure on condensed dust is confirmed. Three equa- 
tions (7) and thus (8), (9) and (10) were obtained in the three 
ranges of effective temperatures distinguished. The usual state- 
ment (van Loon et al. 2003 1 that M oc T~® is approximately true 
only above 2900 K (Eq. (9)) and below 2400 K (Eq. (10)), with 
different proportionality coefficients however. There is practi- 
cally no variation over the 2400-2900 K range in Fig. ^ (in- 
cluding about 55% of our sample). Predictions and observa- 
tions are in agreement provided the variations of luminosity 
on the AGB in the HR diagram (Bergeat et al. \2002b^ and the 
mean masses ranging from about IMq around 3000 K to about 
4Mq around 2400 K (see Table 4), as deduced from theoretical 
tracks in the HR diagram and from pulsations masses (Bergeat 
et al. 2002c% Here mass is the third important parameter, with 
influence much larger than that of the C/O abundance ratio. At 



low effective temperatures, Eq. (10) corresponds to a constant 
mass of about 4Mq. At high temperatures, Eq. (9) describes 
carbon variables of about 1M0 or less, but both observations 
and predictions are uncertain there. Nine giants in our sample 
(at least) do have mass loss rates larger than predicted from 
Eqs. (7) to (10). Four cases of detached shells are found in 
the 2600-2800 K range, i.e. objects having experienced highly 
episodic mass loss possibly caused by a thermal pulse. One of 
them (R Scl) has returned to a moderate mass loss rate (Wong et 
al. 2004 1, and this is certainly true for the others (e.g. Olofsson 
et al. 2000 for TT Cyg). This is a strong argument in favor of a 
specific mass loss mechanism at work. The other five stars are 
very cool objects with T e g < 2500K and optically-thick dust 
shells shown by strong infrared excesses, a consequence of 
very high mass loss rates, like RW LMi (CIT6) or CW Leo 
(IRC+10216). Schoier et al. d2002> have shown that the gas 
mass loss rate of those extreme objects did not vary much over 
the last thousand of years. As argued in Sect. 5, they may be 
stars with masses of about 1.5 - 2.5M0, in contrast with the 
4.2M0 favored for the cool objects with lower rates as given 
by Eq. (10). If true, no specific mass loss mechanism would be 
required here, contrary to the case of detached shells. 

A good agreement is also obtained between predictions and 
observations in the diagram of expansion (outflow) velocity vs. 
effective temperature (Fig.|2j, as described in Sect. 6. The ex- 
pansion (outflow) velocity increases with decreasing effective 
temperature, with about 86 % of the objects within the range of 
predictions as quoted in Table 4 from Eq. (6). To within error 
bars, the observations are located between the two lines ob- 
tained for ec/eo = 1 and ec l^o — 2 respectively, the other pa- 
rameters being kept equal to their values in Table 4 (including 
the mean mass distribution successfully used for mass loss rates 
and described above). This is the range deduced from spec- 
troscopy and model atmospheres (Lambert et al. 1986; see also 
Fig. 7 of Bergeat et al. 2002c 1. In Fig. 12 it is remarkable that 
for effective temperatures higher than about 2900 K, the stars 
concentrate towards the ec/fo = 1 line, below 10-15 km/s. The 
maximum observed velocity rapidly decreases with increasing 
temperature, leaving empty a large triangular area in Fig. [2] The 
result of Bergeat et al. J2001I Fig. 7 in l2002b> is confirmed that 
above 2900 K, the C/O ratio remains close to unity with little 
dispersion. The ec/eo abundance ratio is confirmed here as the 
third major parameter with effective temperature and luminos- 
ity, while stellar mass has little influence in Fig. |2 Cases of 
superwind however blur the information from that figure. We 
conclude that both mass loss rates and outflow velocities are 
correctly modeled from recent hydrodynamical calculations of 
dust-driven winds. Both the levels and the dependence on ef- 
fective temperature are fairly well reproduced. 

Both the predicted and observed ratios of dust-to-gas mass 
loss rates show little variation with effective temperature if 
any. As found in Sect. 7, the average observed value of 
< M d /M g >^ (1.3 ± 1.1) 10~ 3 is a factor 2.2 larger than the 
mean predicted value < pd/g >=< Pd/Pg >— (5-9 + 0.7) 10~ 4 
from Eq. (5). We have discussed in Sect. 7 the influence of the 
stellar mass and C/O abundance ratio. It should be noted that 
the uncertainties are very large here. Among the possible expla- 
nations is a moderate underestimate of abundance ratios (say 
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20%) from model atmospheres, the overal agreement of Fig. [21 
being preserved. Another explanation of the above-mentioned 
disagreement might be some phenomenon in the micro-physics 
of dust grains, especially the non-equilibrium dust formation 
process (Andersen et al. 2003 Sandin & Hofner 2004 and ref- 
erences therein). Allowing drift in the time-dependent (grain 
growth) wind models alters their structure. In some cases there 
is several times more dust in drift models than in non-drift ones, 
even at low drift velocity (a few km/s). It is also found that drift 
plays an active role in accumulating dust to certain narrow re- 
gions (Sandin & Hofner 2004 1, which renders the above com- 
parison even more difficult. It was however seen in Sect. 7 that 
element abundances from model atmospheres raise constraints 
limiting the magnitude of such effects. Further investigations 
are needed before we can determine which estimate is closer to 
truth. 

We finally give in Appendix A a few simple approximate 
formulae in a parametric form. They reproduce the results of 
the present study with fairly good accuracy. They can be use- 
ful for readers wishing not to enter into the full details of the 
present analysis. 
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Appendix A: Approximate formulas 

Some readers may wish not to enter into the full details of the 
present analysis. From Eq. 6 and Fig. |3 we propose for the 
expansion velocities 

logy,,, = (-0.219 + 0.012) log T eft + (1.603 + 0.041) 

+ 1.741og(e c /e ) (A.l) 

for T eff < 3 150 K, and for T eff > 3 150 K, 

log Vco = (-0.682 + 0.056) log T eff + (3.219 + 0.198) 

+ 1.74 log (et/eo) (A.2) 

where the first two terms include the data in Table 4 concerning 
the other (minor) parameters. Concerning the mass loss rates 
(Eq. (4) and Fig.0, we can write 

logM = (-14.62 ± 0.1 8) log T eff + (44.66 + 0.60) 

-2.86 log (M/Mq) (A.3) 

for Teff < 2900, K, and for T eff > 2900 K, 

logM = (-19.12 ± 0.44) log Teft + (60.21 + 1.53) 

-2.86 log (M/M©) . (A.4) 

For the dust-to-gas density ratio, the observations point to 
< Pd/g >=< Pd/Pg >— (1.3 + 1.1) 10~ 3 while predictions from 
Eq. (5) can be written as 

log(p d /Pg) = (1.630 ± 0.268) log Teff + (-9.007 ± 0.906) 

+2.45 log (ec/eo) (A.5) 
for Teff < 2775, K, and for T eff > 2775 K, 
log (p d /p g ) = (4.635 ± 0.605) log Teff + (-19.35 ± 2.10) 

+2.45 log (ec/eo) . (A.6) 



The mean effective temperatures and luminosities (in solar 
units) of Table 4 follow the relation 

logl = (-3.065 + 0.100) log Teff + (14.31 + 0.34) (A.7) 

for Teff < 3050, K, and for T eff > 3050 K, 

log 1 = (-6.528 + 0.693) log T eff + (26.34 + 2.43) (A.8) 

two relations where the true dispersion in the H-R diagram is 
not taken into account. The mass-luminosity relation in Table 4 
was described in detail in Sect. 6 of Bergeat et al. (2002c, their 
equations (10) to (12)). In solar units, we obtain 

logl = (0.72 ±0.07) log m +(3.58+0.12) (A.9) 

for the whole range 0.5 - 4.2 M Q or 3500-2300K. Restricting 
ourselves to the best-documented region from CV2 to CV5 
(3150-2550K) roughly corresponding to 1 - 3 Mq, we derive 

log 1 = (0.497 ±0.075) log m +(3.658 + 0.034) (A. 10) 

that is nearly 

1 oc m 1/2 (A.ll) 
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